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ABSTRACT 

Context. The discovery and chemical analysis of extremely metal-poor stars permit a better understanding of the star formation of 

the first generation of stars and of the Universe emerging from the Big Bang. 

Aims. We report the study of a primordial star situated in the centre of the constellation Leo (SDSS J102915+ 172027). 

Methods. The star, selected from the low resolution-spectrum of the Sloan Digital Sky Survey, was observed at intermediate (with 

X-Shooter at VLT) and at high spectral resolution (with UVES at VLT). The stellar parameters were derived from the photometry. 

The standard spectroscopic analysis based on ID ATLAS models was completed by applying 3D and non-LTE corrections. 

Results. An iron abundance of [Fe/H]=-4.89 makes SDSS J1029I5-I-172927 one of the lowest [Fe/H] stars known. However, the 

absence of measurable C and N enhancements indicates that it has the lowest metallicity, Z < 7.40 x 10"'' (metal-mass fraction), ever 

detected. No oxygen measurement was possible. 

Conclusions. The discovery of SDSS JI029I5-I-172927 highlights that low-mass star formation occurred at metallicities lower than 

previously assumed. Even lower metallicity stars may yet be discovered, with a chemical composition closer to the composition of 

the primordial gas and of the first supemovae. 

Key words. Stars: abundances - Stars: Population 11 - Stars: Population III - Stars: foimation - Galaxy: evolution - Cosmology: 
observations 

1 . Introduction so far. But absence of evidence is not evidence of absence. Such 

O stars, if they exist, could be very rare and faint. Some theories of 

; Very metal-poor stars ai-e the rehc of the primordial Universe. ^^^ formation predict that the stars that formed from the ejecta 

CN _ They are formed from matenal whose composition is close to ^f ^j^^ ^^^^ f^^ generations of primordial massive stars needed 

T: . that of the Universe emerging from the Big Bang, only with ^ chemical pattern different from the one well-known in metal- 

> ; traces of elements heavier than lithium, which were processed p^^^ ^^^.^ ^j^^^^ ^p^^ ^^^^ ^„ enhancement in a-elements, have 

Vp , by the first generations of stars. A low-mass star (M < O.8M0) a scaled solar composit ion (Bonifacio et al. 2009). According to 

^ . formed just after the Big Bang would still shine today on the Bromm & Loeb (S^ and Schn eider et all (12003D, a minimum 

^ ; main sequence, and could be easily recognised from its observed jnetalHcity is necessary to permit the cooling of the collapsing 

- ■ ■ spectrum because only lines of hydrogen and the Li doublet at cloud for small-m ass stars to form. In particular, according to 

670.7 nm would be visible. Such stars have not been detected, |Bromm & Loebl (l20()l . extremely iron-poor stars can form at 

and some models of star formation state that it is impossible .^b-solar masses when enhanced in C and O. The discovery of 

that such star might form. Indeed, it has been speculated that the ^^^^ hyper-iron-poor and very carbon- and nitroge n-enhanced 

first generation of ste^ consisted of massive stars (M > SOM^ ^^^^ in the past decade (iChrisdieb et al. 2002: Frebel et al 



metals that were ejected into space for the following generations (2003'). It is interesting to no te that the coohng bv dust ( Omukai 



IJNon 



r- — ^ — — TT- 1 |, „„ , .„ ' ■ • , sLcus 111 uic pa?ii uc>;au c in^iiiiaiiicu ci ai. z.\j\jz.. ricuci ci ai. 

l Ostnker & Gnedin| [1996) that evolved rapidly and synthetised 2005 UNorris et al.l2007h stipported the theory of iBromm&L^ 



of stars through a supernova type II explosion. This scenario is ^^ ^j 2008; Schneider et al.!|201lD may imply a considerably 



supported by the fact that no zero-metal star has been observed j^^^^ ^^-^-^^^ metallicity, with no specific requirement of chem 
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ical composition. Another fundamental ingredient of star for- 
mation besides cooling is fragmentation. If a large-mass cloud 



Based on observations obtained at ESO Paranal Observatory, GTO ^ ^ • ^ n • ^u ..1. ■ r 

„„, „ „„„. , oo^ T^ crv^c can fragment into smaller pieces, then these pieces can form 

programme 086.D-0094 and programme 286.D-5045 , ° . , K r , , ■ 

j^* „,. T- 11 low-mass stars even in the absence of metals, the main cooling 

Gliese Fellow '^ 
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agent beir ig molecular hydrog en (iNakamura & Umemural200l]) . 
Recently, IClark et aP (1201 l l) have shown that the first genera- 
tion of stars are indeed likely to form as binary or multiple sys- 
tems resulting from the fragm entation of an initia l cloud of large 
mass. Under these conditions. iGreif et al.l (1201 ll) derived an ini- 
tial mass function (IMF) that is essentially flat between 0. 1 and 
10 solar masses. This IMF is top-heavy, in the sense that most of 
the mass is concentrated in a few massive stars, and yet it allows 
for the formation of low-mass stars. In this latter scenario there 
is no critical metallicity. 

We present here the analysis of SDSS J102915-H172927, an 
extremely metal-poor star, that shows no strong enhancement in 
carbon and nitrogen, and has a global metallicity more than four 
orders of magnitude lower than the solar one: Z < 7.40 x 10^^, 
to be compared to Zq = 1.53 x 10"^ (ICaff'au etal.ll2011ah . giving 
a metallicity in respect with solar one of Z = 5 lO'^Z©. If CNO 
abundances are non-enhanced, Z - 7.15 x 10~^. With respect t o 
the results presented earlier by our group (ICaffau et al.ll201 la) , 
the present analysis uses of a custom CO^BOLD hydrodynami- 
cal model and NLTE computations for Mg, Si, Ca, Fe, and Sr. 



2. Observations and data reduction 

For SPSS J 1029 15 -Hi 7 2927 we av ailed ourselves of X -Shooter 
(lD'Odoricoetal.ll2 006') and UVES (iDekke r et al.'2000') observa- 
tions. The star was observed during the X-Shooter French-Italian 
Guaranteed Time of Observation (GTO) on 10 February 2011, 
programme 086.D-0094, P.I. P. Bonifacio. The spectral range 
covered is 330-2400 nm, split into three spectral arms (UVB, 
VIS, and NIR). For the integral field unit (IFU), which re-images 
an input field of 4"x 1.8" into a pseudo slit of 12"x 0.6" and a 
1x2 binning, staring mode was used. Using the IFU as a sheer, 
the resolving power was R=7900 in the UVB arm and R= 12600 
in the VI S arm. The spectra were reduced using the X-Shooter 
pipeline dGoldoni et al ."2006). The details o n X-Shooter obser- 
vations can be found in iCaffau et al] (1201 Ibl) . 

On the basis of the highly promising analysis of the X- 
Shooter spectrum we applied for ESO Director Discretionary 
Time (DDT) to observe the star with UVES@VLT, and the ob- 
serving blocks (OB) were observed on March- April 201 1, within 
the programme 286.D-5045, P.I. P. Bonifacio. Five OBs were ex- 
ecuted in Die #1, standard 390nm + 580nm setting, two in DIC 
#2, standard 473nm + 760nm setting. All exposures had a length 
of 3000s, slit width was set to 1 .4", which led to an expected res- 
olution of ~38 000. However, during most exposures the seeing 
was better than 1 .4", thus leading to slightly better resolution in 
the final spectra. The signal-to-noise ratio per pixel at 650 nm 
was between 30 and 40 for all spectra. 

3. IVIodel atmospheres 

The analysis was mainly performed with ID plane-par allel hy- 
dros tatic model atmospheres, computed with ATLAS 9 (' Kurucj 
|1993,.2005.) in its Linux version (.Sbordone.2005.) . The line opac- 
ity was treated through opacity dist ribution functions (ODFs) , 
and we used the ODFs computed by ICastelh & Kuruczl ( l2003l) 
for a scaled solar metallicity corresponding to [Fe/H]=-4.5. The 
ATLAS 9 model was computed with the stellar parameters de- 
rived from photometry, as given in Tab.[T] The 3D abundance 
corrections were derived with a 3D-C0^B0LD model expressly 
computed for this star and the related IDlhd model, computed 
with the LHP code (for details concerning the 3D corrections see 
ICaffau et al.ll201 la) . To derive the 3D corrections on the chem- 
ical abundances we compared the abundance derived by using 



Table 1. Characteristics of the star SDSS J102915H-172927. 



Star 


SDSS J102915+ 172927 


RA (J2000.0) 


10/! 29/71 15. 15 J 


Dec (J2000.0) 


-l-17°29'28" 


/ 


221.785° 


b 


-1-55.863° 


Epoch of SDSS photometi-y 


2005.929 


u, SDSS 


17.736 ±0.011 


g, SDSS 


16.922 ± 0.004 


r, SDSS 


16.542 ± 0.004 


i, SDSS 


16.388 ±0.004 


z, SDSS 


16.330 ±0.008 


J, 2MASS 


15.513 ±0.052 


H, 2MASS 


15.126 ±0.077 


K, 2MASS 


15.146 ±0.116 


Av 


0.084 


V,,d, [kms-'] 


-34.5 ± 1.0 


Teff, [K] 


5811 ±150 


log 9 


4.0 ± 0.5 


^[kms-'] 


1.5 


Metallicity, SDSS 


-3.73 


[Fe/H] 


-4.89 ±0.10 


[«/H] 


0.23 ± 0.26 



the 3D-C0''B0LD model with the abundance derived with ID 
plane-parallel model, the IDlhd model. The IDlhd model em- 
ploys the same micro-physics and radiative transfer scheme as 
the 3D-C0^B0LD model, so that any difference in the abun- 
dance determination is due to the 3D versus ID treatment of 
convection. The 3D model used in the following was com- 
puted for stellar parameters Teff - 5850 K, loggf= 4.0, [Fe/H] 
- -4.0. The computational domain is a cartesian box of size 
25.8 X 25.8 X 12.5 Mm^ resolved by 200 x 200 x 200 ceUs. 
The frequency dependence of the total radiative opacity was 
approximately taken in to account by adoptin g the opacity bin- 
ning method (see e.g. 'Frevtag et al.l l2012h . using a total of 
1 1 frequency groups. The detailed monochromatic (continuous 
plus line) opacities needed to evaluate the average opacity for 
each frequency group as a function of P an d T were derived 
from the M ARCS stellar atmosphere package (iGustafsson et al.l 
I2OO3LI2OO81) formetaUicity [Fe/H]=-4, [a/H]=0.4 (B. Plez, priv. 
comm.). 

The present 3D model differs from the one pre viously used 
for the calculation of 3D abundance corrections in ICaffau et al.l 
(1201 id) in several aspects: (i) its metallicity is lower by a factor 
of 10 ([Fe/H]=-4 instead of -3), (ii) it uses a better frequency 
resolution (11 instead of 6 frequency groups), (iii) it has a higher 
spatial resolution (200 x 200 x 200 instead of 140 x 140 x 150 
grid cells). Nevertheless, the mean temperature stratifications 
are quite similar (see Fig.[T}, except for the upper photosphere 
(log tross < -3) where - contrary to expectations - the lower 
metallicity model shows somewhat higher temperatures. Both 
3D models agree in the amplitude of the photospheric tempera- 
ture fluctuations, which is fairly low, presumably because of H2 
molecule formation in the higher photosphere. The most striking 
difference between the 3D and the ID models is the temperature 
in the outer photosphere, which in both 3D models is signifi- 
cantly cooler than their IDlhd counterparts. 

The total 3D LTE abundance corrections (A(3D) - 
A( IDlhd)) were computed for all atomic lines measured in the 
observed spectra, assuming q'mlt=0.5, and a microturbulence 
parameter^ = 1.5 km s"' for the IDlhd models. For the molec- 
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Optical depth 



Fig. 1. Mean temperature structure as a function of Rosseland 
optical depth of the CO^BOLD model used in this work (Teff 
=5850 K, log 0=4.0, [ Fe/H1=-4 .0; thick [blue] line), and of the 
3D model used in Caffau et alJ ^OMlc) (5850/4.0/ - 3.0; thin 
[red] line). For each model, the mean temperature is obtained by 
averaging the temperature as yiT^ over iso-TRoss surfaces of 20 
selected snapshots. Shaded areas bounded by dashed lines indi- 
cate the ±lcr range of deviations. Small open symbols outline 
the temperature stratification of the IDlhd model with [Fe/H]=- 
4. 



ular bands, G-band and NH-band, representative lines were con- 
sidered. Obviously, the physical reason for the large negative 3D 
corrections for Fe i and Ni i is the large 3D- ID temperature dif- 
ference in the upper photosphere. 



4. Analysis 

The reference solar abundances we used were taken from Caffau 



et al. (1201 la!) for C, N. O. and Fe. and from lLodders et al.l(l2009h 
for the other elements. The atomic data for the lines were taken 
from the large programme "First Stars" lead by Roger Cayre l 
dBonifacio et alJl2009l : ICavrel et al.ll2004l:lFran9ois et al.ll2007l) . 

4.1. Stellar parameters 

We der ived the effective te mperature of 58 11 K from the (g - z) 
colour (Ludwig et al. 2008). We also fitted the Hff wings. By us- 
ing a grid of synthetic spectra computed from an ATLAS model 
and a modified version of th e BALMER codq3 that uses the the- 
ory of Barklem et al.' ('2000'.^ for the self-broadening and the 
profiles of Stehle & Hutcheon (1999) for Stark broadening, we 
derived an effective temperature within few K. However, when 
using a grid of 3D profile the temperature is about lOOK hotter. 
We assumed an uncertainty of 150 K in the effective tempera- 
ture. This value is derived by linearly adding the uncertainties in 
the SDSS (g - z) colour and an estimate on the error on the zero 
point of the calibration. We took 75 K as the uncertainty due to 
the eiTor in the colour on the derived Teu , which corresponds to 
three times cr=25 K. Comparing the r eff derived from the (a - z) 
colour of the four calibrator stars of Fukugita et al.l (1 19961) and 
the IRFM temperatures of Alonso et al. (1996), we find an av- 
erage difference of about 70 K. A similar comparison with the 



' The original version 
|http : //kurucz . harvard . edu/i 



IRFM temperatures of lGonzalez Hernandez & Bonifaciol (l2009h 
provides a difference of the order of 10 K. Conservatively, we 
adopted the difference with the Alonso et al. (1996) tempera- 
tures as an estimate of the zero point error of the calibration. 

UVES spectra were able to give us better insight than the 
X-Shooter data in the gravity determination. The Ca ionisation 
equilibrium, when computed in LTE, supports the adopted grav- 
ity (loggr= 4.0). By comparing the (m - g) colour with the the- 
oretical colours computed from the grid of ATLAS 9 models of 
ICastelU & Kurucz I j2003lFl we derived loggf= 3.9. The use of 
the u - g colour for gravity determination implies a fairly large 
uncertainty, both for the sensitivity to photometric errors (an er- 
ror of 0.05 mag in the colour translates into an error of 0.3 dex 
in log g) and for the uncertainty in the zero point. A comparison 
of the u - g derived gravities with the Hipparcos parallax based 
gravities of Gonzalez Hernand ez & Bonifacio (200 9) for the four 
SDSS photometry calibrators (iFukugita et al.ll996h suggests that 
the u - g colour provides gravities that are systematically lower 
by 0.5 dex. We cannot exclude that the star has indeed a higher 
gravity of about 4.5. The non-LTE Cai/Caii equilibrium would 
require an even higher gravity of about log g= 4,8, which is diffi- 
cult to accept considering the gravities derived from isochrones, 
which are compatible with the effective temperature of the star. 
Nevertheless, all this evidence allows us to robustly exclude that 
the star is either a sub-giant or a horizontal-branch star. 

No Fe II line is detectable in the spectrum, therefore the iron 
ionisation equilibrium cannot be adopted for gravity determina- 
tion. But lines of both ionisation states of Ca are present in the 
UVES spectrum. The ID-LTE abundance derived from the three 
lines of the IR Ca ii triplet differs by 0. 1 2 dex from the abundance 
derived from the Cai line at 422.6 nm in the ID-LTE analysis. 
This latter line is insensitive to gravity, while the IR Ca ii lines 
depends on the choice of gravity. This good agreement between 
the Ca lines could be an indication that the gravity is correct, 
but it is broken when we apply non-LTE corrections. In this case 
the difference becomes 0.38 dex. The use of Sh = 1 instead of 
the adopted value of 0. 1 (see section l4~8] ) cannot remove the dif- 
ference between the two ionisation stages, although it does de- 
crease (0.23 dex). The higher Ca abundance from the neutral Ca i 
line would require higher gravity. When considering a gravity of 
loggf=4.5, the Ca abundance from the Cai line remains mostly 
the same (within 0.03 dex) both in ID-LTE and ID-NLTE, while 
the abundance of Ca derived from the IR-triplet would increase 
by 0. 12 and 0. 17 dex in ID-LTE and ID-NLTE, respectively. The 
agreement would definitely improve, but the increase in gravity 
would not be enough to result in an agreement of the Ca abun- 
dance derived from the two ionisation states of Ca, and we can- 
not support a higher value of gravity in this star. To determine 
the gravity we also fitted the wings of high-order Balmer lines 
(He and higher) and found that the gravity should be in the range 
4.3 <log g< 3.6 when the temperature of the synthetic spectra is 
fixed at 581 1 K. By increasing the temperature by about 200 K, 
the gravity increases by 0.5 dex. From all these results we de- 
cided to keep log gf=4.0 as our value for the gravity, but the star 
could be a main-sequence star with log ^=4.5 or a turn off-star 
with logg> 3.5. We are confident that the star is no sub-giant. 
We cannot exclude that the problem with the Ca abundance is 
related to the effective temperature. An increase of +200 K in 
effective temperature would change the calcium abundance de- 
rived from the 422.6 nm line by about ±0.2 dex, and the abun- 
dance derived from the Ca ii lines of the triplet would change by 
about +0. 1 dex. 
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Fig. 2. Best solutions of the fitted isochrones to the eight-band 
SDSS-2MASS photometry of SDSS J102915+172927. In the 
upper panel we depict the solution residuals. 



We deri ved the micro-turbu lence of 1.5 km s' by using the 
formula in lEdvardsson et al.l ([1993). The same value of the 
micro-turbulence was applied to the IDlhd model used for com- 
puting the 3D corrections. The 3D corrections are sensitive to the 
choice of the micro-turbulence, because the abundance from the 
IDlhd model depends on it in the same way as the abundance 
derived from the ATLAS model. The ID abundance corrected by 
the 3D corrections is therefore largely insensitive to the choice 
of micro-turbulence. 

The radial velocity of the star from the UVES spectrum is 
-34.5 km s \ which agrees well with the radial velocity derived 
from X-Shooter and SDSS spectra. 



4.2. Variability 

We investigated possible long term-photometric variabilities of 
SDSS J102915+172927 and, given the material at our disposal, 
we can exclude this at a level of Ar < 0.006. A detailed descrip- 
tion can be found in the appendix. 



4.3. Distance (isociirone fitting) 

We fitted the eight-band photometry (SDSS-I-2MASS) as re- 
ported in Table[I]of SDSS J102915+172927 with different sets 
of low-metallicity isochrones. There are no published isochrones 
that match the low metallicity of this type of object. We col- 
lected a sample of the most re cent low-metalli city isochrones 
and added two special cases: Mar igo et al.l (l2003h who presented 
isochrones for zero-metallicity low-mass stars but limited to 
0.70Mq; and an unpublished set of FR ANEC isochrones, com - 
puted according to the prescriptions of IStraniero et al.l (Il997h . 
kindly provided by Chieffi and Limongi, with [Fe/H]= -6.0. The 
sets of all theoretical isochrones used are listed in Table |2] 

The fitting was made using all photometric measurements 
from u to K weighting on the photometric precision. This is 
particularly important since the infrared 2MASS colours of this 
star are at the lower limits of detection and should be treated 
with care. For each stellar mass of the isochrone we varied only 
the distance modulus to calculate x^- We then chose the best 
model on the basis of the lower ;^'^. Considering that the grav- 



ity from spectroscopy is compatible with a dwarf star, we as- 
sumed that the star is on the main sequence. The age range as- 
sumed is 12.0 or 13.0 Gyr depending on the availability of the 
isochrones. We checked the goodness of the solutions also along 
the sub-giant branch, but we always found that the solutions 
were worse than the dwarf solutions. All isochrones used are 
on the VEGA Johnson-Cousins system, except for the Girardi 
isochrones (Girardi et al. 2002, 2004, 2005), which can be ob- 
tained for a variety of different photometric systems and for 
which we chose SDSS+2MASS directly. To obtain all magni- 
tudes on this system, we properly transformed the isochrones 
in the observed system using the colour transformations pro- 
vided by the SDS!^ This of course has some effect on the given 
precision of the solutions. In the process of fitting the photom- 
etry, each isochrone was reddened with the observed value of 
extinction given by the Schlegel maps (see Table [TJ using the 
proper extinction coefficients for the SDSS/2MASS passbands. 
The output of the whole procedure is the value of the distance 
modulus and the stellar mass of the isochrone. 

The calculated best solution together with the residuals are 
plotted in Fig.|2]where we show the five isochrones used in this 
study. In Table |2] we list the parameters for each of the best so- 
lutions together with the x^ value. In the figure the solutions 
overimposed on the photometric measurements are shown in the 
lower panel, while the residuals of the fitting are plotted in the 
upper panel. While for four isochrones we obtained the abso- 
lute best solu tions for th e distan ce and mass, for the Maiigo et 
al. (2003) data our search stopped at the minimum stellar mass 
available, 0.70 M©, but clearly the fitting trend required a lower 
stellar mass of the order of ^ 0.60 Mq. Future versions of these 
isochrones (Bressan et al., in preparation) will allow a more pre- 
cise constraint of the evolutionary status of this star. 

Finally, it is not surprising to find that the best isochrone is 
that of Girardi et al. because it was "coloured" directly in the 
SDSS/2MASS systems. Our preference is the Chieffi& Limongi 
isochrone because it gives the lowest ;^f^, the best match with 
the spectroscopic metallicity, and with the Teff of the star (a dif- 
ference of only ^ 35° K). To summarise, the assumed distance 
modulus is (m - M)o = 10.68 + 0.15 mag or 1.37 + 0.20 kpc. 

4.4. Kinematics 

Kinematical data are fundamental for understanding the nature 
of SDSS J102915-H172927. Radial velocity and proper motion 
measurement allowed us to compute a Galactic orbit. 

4.4.1. Radial velocity 

The radial velocity of the star as measured from the summed 
spectrum of all UVES data is -34.5 + 1.0 kms"', which agrees 
well with the -31.0 ± 10.0 kms^'measured from the X-Shooter 
spectrum. We found no significant variations of the radial veloc- 
ity. The difference of UVES and X-Shooter is well within the 
given error estimates. Given the low value measured, we per- 
formed an analysis of the compatibility of the radial velocity 
with the different stellar populations of the Galaxy. We used the 
Besangon Galactic model ( Robin et al. 200 3) to simulate a field 
of 1 square degree surrounding the star. In the simulated field we 
selected a sample of stars with similar magnitude and colours 
around SDSS J102915-I-172927: we used a box 0.1 mag wide in 
both r and (r - /)■ The observed radial velocity lies within 1.2cr 
from the expected mean radial velocity of the field stellar pop- 



http://www.sdss3.org/dr8/algorithms/sdssUBVRITransform.php 
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Table 2. Parameters of the best-fitted (x^) isochrones for main-sequence stars. 



Author 


[Fe/H] 


Y 


[tt/Fe] 


Age 


Mass 


logs 


Teff 


(m-M)o 


X 








[dex] 


[Gyr] 


[Mo] 




[K] 


[mag] 




Marigo et al. 


Zero 


0.249 


_ 


12.6 


0.70 


4.6 


6108.0 


11.09 


0.070 


Chiefli & Limongi 


-6.00 


0.232 


- 


12.0 


0.68 


4.7 


5775.3 


10.68 


0.029 


Pietrinfemi et al 


-3.26 


0.245 


+0.40 


12.5 


0.64 


4.6 


5751.5 


10.51 


0.033 


Yale-Yonsey 


-3.76 


0.230 


+0.60 


13.0 


0.62 


4.7 


5592.4 


10.29 


0.038 


Girardi et al. 


-2.45 


0.245 


+0.00 


13.0 


0.67 


4.6 


5873.5 


10.76 


0.011 
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Fig. 3. Measurement of the proper motion in the two coordinates, 
RA in the bottom panel, DEC in the upper panel. The zero points 
for the offsets are given by the position in the SDSS-DR7 cata- 
logue. The results of the two different fits are shown: for all po- 
sitions (long dashed line) and for the positions since 1990 only 
(dotted line). 



ulation, thin and thick disc and halo, which has an average of 
+37.4 kms"' and a cr = 58.2 kms"'. A "normal" radial veloc- 
ity excludes the possibility that this star is a member of recently 
accreted streams or other substructures present in the halo of the 
Galaxy. 



4.4.2. Proper motion 

Although recent surveys and old photographic catalogues have 
allowed us to build all-sky catalogues of proper motions for sev- 
eral millions of stars, subtle systematic errors are always lurking 
in the data, therefore we recommend caution. 

In the present case the r ecent PPMXL Catalogue (" Roeser et 
al. |201(|) provides a proper motion measurement for this star: 
fj.aCos(6) = -12.8 + 3.9 mas/yr and fis = -6.7 + 3.9 mas/yr. 
The values are based on six different measurements comprising 
POSS I and II, 2MASS, and SDSS. This is not the only available 
measurement since also SDSS-DR7 Munn et al. (2004) provides 
a measurement for the proper motion: yL(„cos((5) = -19+2 mas/yr 
and yu5 = -1+2 mas/yr There also seems to be a negative value 
for the right ascension here, but almost no motion in declination. 
However, the latter measurements were discovered to have sys- 
tematic errors (Munn et al. 2008) that were coiTected, but leave 
some doubt about the derived proper motions. We have to note 
that the faintness of the object makes it intrinsically difficult to 



measure in some of the surveys (like in 2MASS and UCAC3) 
because it is near the magnitude limit of the detectors. 

To check the validity of the above values of proper mo- 
tions, we independently measured the positions of the star in 
images of different epochs. We then combined these positions 
with the star's positions in different recent astrometric catalogues 
based on CCD detectors. The proper motion strongly depends 
on the first epoch position given by the POSS-I plates taken in 
the 1950s. Since then there is a long gap until the 1980s. To 
improve the first epoch position, we recalculated the astrome- 
try of the available photographic plates, POSS-I and POSS-II 
with the help of the SDSS photometry in the field. We recovered 
1x1 deg- images for both POSS-I/II around our target star in 
the blue and red from the MAST database at STScI. We then 
obtained a new astrometric solution for each POSS image us- 
ing the SCAMP programmeQ The extracted sources were com- 
pared with an SDSS-DR7 reference catalogue built by select- 
ing only objects classified as compact galaxies with magnitude 
r < 21.00. A total of 2532 galaxies were used as reference out 
of a total of 20532 objects present in the field. With this sample 
of galaxies we minimised the systematic effects present in in- 
dependent astrometiies. We then "swaiped" each image with the 
SWARP programme (see footnote|4]i and measured the positions 
with SEXTRACTOR (parameters X/Y-WIN JMAGE and given 
errors). We found no significant presence of systematic trends 
either with position in the plate, magnitude, or type of object. 
We verified that the used reference galaxies defined a common 
zero point of the two astrometries (POSS and SDSS-DR7). 

A similar procedure has been used also on the recent R- 
and V-band images taken at the Asiago Cima Ekar Schmidt 
Telescope already described above. We used the same refer- 
ence list of galaxies to perform the astrometrisation of the field, 
again using the SCAMP/SWARP/SEXTRACTOR sequence of 
programmes. In this case we started for each band from the four 
images and "swarped" the resulting final positions with sextrac- 
tor to a single image. 

To the 2 POSS-I, the 2 POSS-II, and two Asiago positions 
we added 2MASS, UCAC3, CMC14 and, of course, SDSS-DR7 
to compute the final proper motions that are shown in Figure |3] 
Clearly, the entire "signal" of the proper motion is indeed given 
by the first epoch positions of the 2 POSS-I plates. Fitting linear 
regression gives a solution of fiaCos{S) - -12.33 + 0.66 mas/yr 
and fig = -6.60 + 0.99 mas/yr, which is quite similar to the 
PPMXL solution. If we drop the POSS-I and use only the mea- 
surements taken after 1990.0, we find a less significant solution 
fia cos(<5) - -7.23 + 1 .96 mas/yr and fig - -5.01 + 3.00 mas/yr. 



E. Bertin, see for details http://www.astromatic.net/ 
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STAR J 102915,1 4 





3.0 



Fig. 4. Galactic orbit for star J 1 029 1 5 . 1 4 + 1 72927 . 9 obtained us- 
ing the PPMXL proper motion. The top-left panel shows the X/Z 
plane; top-right shows the Y/Z plane; bottom-left shows the X/Y 
plane; bottom-right shows the meridional plane R/Z. 



4.4.3. Galactic orbit 

Combining the positional and kinematic information, we calcu- 
lated the Galactic orbit of the star by performing several tests 
varying the input parameters within the errors. The code used 
(Gratton et al. 2003) is based on the Allen & Santillan (1991) 
prescriptions for the Galactic potential components with updated 
parameters for the constants describing the potential. Although 
the low height on the Galactic plane (Z ~ 1.0 kpc) may suggest 
a thick disk orbit, this can be safely ruled out. The orbit solution 
indicates that the star belongs to the halo with the maximum 
height above the galactic plane Z^ax = 4.8 + 0.4 kpc, the orbital 
apocenter at Rmax = 9.6 + 0.6 kpc, and is plunging towards the 
Galactic centre, with orbital pericenter R^m = 0.9 + 0.1 kpc. See 
FigSl Adopting the proper motion values obtained in the previ- 
ous section from the positions after 1990.0, we obtain a similar 
orbit with a more extreme orbital pericenter R,nin - 0.4+0.1 kpc. 
An even more extreme value of 0.2 ±0.1 kpc is obtained when 
we adopt a null value of the proper motion. 

4.5. Abundances 

Very few lines are measurable in the X-Shooter spectrum. The 
Mg i-b triplet is not visible. Of the IR Ca ii triplet lines, only 
the one at 854.2 nm is clearly visible, but is contaminated by a 
feature produced by the sky subtraction. Some Fe i lines can be 
guessed, not really measured. The only clearly detectable line is 
the Ca ii-K line at 393.3 nm. Its EWof49.2pmis consistent with 
an abundance of [Ca/H]=-3.9. But the measured radial velocity 
is -30km s', comparable to the X-Shooter UBV arm resolu- 
tion of 7 900, which suggests that the line is contaminated by 
the component from the interstellar medium (ISM). From the 
X-Shooter spectrum, we can deduce that this spectrum belongs 



2.0 



1.0 



0.0 




Call-H 



Call-K 





388 



390 



392 394 

Wavelength [nm] 



396 



398 



Fig. 5. Range of the Can H and K lines. From top to bottom, 
SDSS, X-Shooter, and UVES spectrum (solid black), overim- 
posed the synthetic profile with metallicity -4.5, a-enhanced by 
0.4 dex (solid green). 



to an extremely metal-poor star and put an upper limit on the 
metallicity of about -4.0 with respect to the solar metallicity. 

The UVES spectrum resolves the stellar and IS components 
of the Can-K and Caii-H line (see Fig.|5]l. The EW of the stel- 
lar Ca ii-K line is of 27.7 pm, corresponding to an abundance of 
[Ca/H]=^.47. We did not take this line as an abundance indi- 
cator, because it is difficult to separate the stellar and IS compo- 
nents, and because with an EW of 27.7 pm this line is saturated 
and not very sensitive to abundance. 

In the UVES spectrum we can see lines of iron peak elements 
(Fei, Nil) and a-elements (Mgi, Sii, Cai, Can, Tin). For the 
light elements Li and C-N we were unable to find an evident 
signature in the spectra, therefore we can provide only an upper- 
limit. 

For the abundance determination we relied on line profile 
fitting because some lines happen to be blended (sometimes sev- 
eral lines of the same element) and some lines lie on the wings of 
hydrogen lines. We computed a grid of synthetic spectra with the 
effective temperature and gravity of the star, varying in [Fe/H] by 
0.2 dex. We fitted the Fe i features to derive the ID-LTE [Fe/H]. 
To derive the abundances of the other elements, we computed 
grids of synthetic spectra with fixed [Fe/H], by varying the abun- 
dance [X/Fe] of the element X by 0.2 dex and then by fitting 
the line profiles. The abundances derived for all elements are 
repo rted in Table|3] The differences with respect to Table 1 of 
Ca ffau et al.l (l20IId) are that we use a different 3D model atmo- 
sphere and the applications of NLTE. 
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4.6. The Li abundance 



A 3D-NLTE JSbordone et al.ll2010l) Li abundance of 2.2 (Spite 
plateau) would imply an EW for the Li doublet at 670.7 nm of 
about 4.7 pm in this star. Such a feature should be visible in the 
observed spectra, but no sign of the line is detectable in the ex- 
pected wavelength range. In the X-Shooter spectrum, taking into 
account its signal-to-noise ratio ( S/N) and the r esolution, we ex- 
pect according Cayrel's formula (ICavrellll988l) that the limit for 
a feature to be visible is about 1 .5 pm (3(t), which would corre- 
spond to a A(Li)=1.7, close to the Li abundance derived for the 
cooler component of the bin ary system CS 22876-32 ( Gonzalez 
Hernandez et al. 2008). From the S/N of the UVES spectrum 
(160) an upper limit on the EW of 0.1 pm implies A(Li)< 1.1 at 
5cr or A(Li)< 0.9 at 3o-. 

This implies that the star is far below the Spite plateau. This 
may be linked to the fact that at extremely low metallicities 
the Spite plateau displays a "meltdown" (Sbordone et al. 2010), 
i.e. an increased scatter and a lower mean Li abundance. This 
meltdown is clearly visible in the two components of the ex- 
tremely metal-poor binary system CS 22876-32 ([Fe/H]=-3.6, 
the primary with an effective temperature of 6500 K, the sec- 
ondary of 5900 K), wh ich show a different Li c ontent ( Gonzalez 
Hernandez et al. 2008D- The primary lies on the Spite plateau, 
while the secondary lies below at A(Li)= 1.8. The reasons for 
this meltdown are not understood. It has been suggested that a Li 
depletion mechanism, whose efficiency is metallicity-dependent, 
might be able to explain the observations. If this were the case, 
the Li abundance in SDSS J 1029 15 -Hi 72927 would result from 
an efficient Li depletion caused by a combination of extremely 
low metallicity and relativ ely low tem perature. If the star were 
a horizontal branch star (IHansen et al.ll20III) . it would be nor- 
mal for it to be Li depleted. However, we already argued that 
low gravities, compatible with an HB status, are ruled out. A 
sub-giant status should not imply a high Li depletion. The ab- 
sence of Li could be expla ined if SPSS JI02 9 15+ 172927 were 
a "blue straggler to be" (iRvan et al.ll2002h . In this case we 
would expect a measurable line broadening owing to rotation. 
In our UVES spectra we cannot derive any line broadening 
above what is caused by the instrumental resolution, which is 
set by the seeing. Therefore all available evidence suggests that 
SDSSJ102915+172927 is in an evolutionary status from the 
main sequence to the sub-giant branch. 

4. 7. CNO upper limits 

No strong carbon enhancement is evident in the UVES spectrum 
(see Fig.|6]l. We fitted the range of the G-band at 430 nm, and 
found [C/H] ID = -3 .8 1 . The S/N in the range is of about 70. The 
same fit of a synthetic spectrum with S/N of 50 gives [C/H]=- 
4.25, C-enhancement of 0.25 dex with respect to the input value. 
We consider [C/H]< -3.81 as an upper limit. If we applied 3D 
coiTections, the abundance decreased and the upper limit be- 
came more strict. This is because the temperature structure of 
a 3D model is cooler than the reference ID model in the ex- 
ternal layers. In these cooler layers, the formation of molecules 
is favoured, which means that with the same abundance the ob- 
served molecul ar bands are deeper which in turn implies a lower 
ab undance. seelGonzalez Hernandez et al.l (12010) and Behara et 
al. (12010D for details. 

For the nitrogen we investigated the NH band at 336 nm. 
We found no evidence of enhancement (see Fig.|2l) here either; 
when we fited the range we obtained an upper limit of an N- 
enhancement of 0.24 dex. As for the G-band, we consider this 



HE 1327-2326 



SDSS J102915+172927 





430 432 

Wavelength [nm] 



Fig. 6. G-band range of the the star SDSS J 1 029 1 5 H- 172927: the 
observed spectrum (solid black) is compared to the known C- 
enhanced star HE 1327-2326. A zoom on the G-band is shown 
(solid black) compared to the best fit of the range (solid red) and 
with a synthetic profile enhanced in C by 2 dex (solid green). 
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Fig. 7. NH-band range of the the star SDSS JI02915H-172927: 
the observed spectrum (solid black) is compared to the best fit of 
the range (solid green) and with a synthetic profile enhanced in 
N by 2 dex (solid green). For comparison the known N-enhanced 
star HE 1327-2326 is shown. 



value as an upper limit for nitrogen. The 3D corrections would 
imply a lower upper limit for this molecule as well. 

SDSS J102915+172927, with a metallicity of Z = 5 10"%, 
is the presently known object whose composition closest resem- 
bles the primordial one. We computation the metallicity with the 
ID abundances derived from the UVES spectrum, the IDC and 
N upper limits that we derived by fitting the molecular bands 
of CH and NH, respectively, and assuming an enhancement in 
oxygen of [O/Fe]=+0.6. 

4.8. Non-LTE effects on fhe element abundance 
determinations 

Our present investigation is based on the non-local thermody- 
namic equilibrium (non-LTE) line formation for six chemical 
species. We used the non-LTE methods treated in our earlier 
studies and documented in a number of papers, in which atomic 
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Table 3. SDSS J102915+172927. The abundances [X/H] from the line profile fits are given. The adopted log gf values are given 
in Tabl^HU 



Element 






XAillD 




N lines 


Sh 


A(X)o 




ID-LTE 


+3Dcor 


+NLTE 


+ 3Dcor + NLTE 








C 


<-3.8 


<-4.5 






G-band 




8.50 


N 


<-4.1 


<-5.0 






NH-band 




7.86 


Mgi 


-4.71 ±0.11 


-4.68 ±0.11 


-4.52 ±0.11 


-4.49 ±0.12 


5 


0.1 


7.54 


Si I 


-4.27 


-4.30 


-3.93 


-3.96 


1 


0.1 


7.52 


Cai 


-4.72 


-4.82 


-4.44 


-4.54 


1 


0.1 


6.33 


Can 


-4.60 ±0.11 


-4.73 ± 0.03 


-4.82 ± 0.02 


-4.95 ± 0.09 


3 


0.1 


6.33 


Tin 


-4.75 ±0.18 


-4.83 ±0.16 


-4.76 ±0.18 


-4.84 ±0.16 


6 


1.0 


4.90 


Fei 


-4.73 ±0.13 


-5.02 ±0.10 


-4.60 ±0.13 


-4.89 ±0.10 


43 


1.0 


7.52 


Nil 


-4.55 ±0.14 


-4.90 ±0.11 






10 




6.23 


Srn 


<-5.10 


< -5.25 


< -4.94 


< -5.09 


1 


0.01 


2.92 



data and the problems of line formati on were consider ed in 
detail, i.e., in Izhao et al.l (TlQQS. Mgi), IShi et al.l (I2008L Siil 
iMashonkina et alJ (l20oTcai-Caii) Mashonkina et al.1 (1201 iL 
Fei), and iBelvakova & Mashonkinal (il997 . Srii). To solve the 
coupled radiative transfer and statistical equilibrium (SE) equa- 
tions, we used a revised version of the DETAIL programme 
dButler & Giddin^s1ll985h based on the accelerated lambda iter- 
ation, whic h follows the efficient method described bv Rvbicki 
& Hummer (l l99lLll992 '). The non-LTE calculations were per- 
formed for the ATLAS 9 atmospheric structure. 

It is worth noting that we used the most accurate atomic 
data for each chemical species. Photoionisation cross-sections 
are f rom the opacity project (OP) calculations (ISeaton et al.l 
11994 ). Their accuracy is estimated at the level of 10%. This 
small uncertainty translates into an abundance error of no more 
than 0.01 dex. When no OP data are available for Sr ii, we ignore 
photoionisation since it affects its SE only weakly because Sr in 
is an extremely small fraction of the total strontium atoms. The 
atomic data for Ca ii were updated by applyin g effective colh- 
sion st rengths from the 7?-matrix calculations of lMelendez et al.l 
(|2007|) . We accounted for inelastic collisions both with electrons 
and neutral H particles in SE calculations. Hydrogen collisions 
were computed using the formula of Steenbock & Holwegeg 
d 19841) with a scaling factor S-a— 0.1 and 1. For Cai-Caii and 
Fe i-Fe ii our favourite is S'h = 0. 1 as estimated empirically from 
the different influence of hydrogen atom collisions on the differ- 
ent lines of a given ato m in solar a nd stella r spectra ( Mashonkina 
et al. 2007. 2011) . For Srii, we applied ^h = (no hydrogen ic 
collisions), as recommended bv lMashonkina & Gehrenl(l2001h . 

The main non-LTE mechanism for the minority species in 
the model 5811/4.0/-4.5, Mgi, Sii, Cai, and Fei, is the overi- 
onisation caused by superthermal radiation of non-local origin 
below the thresholds of the levels with fiexc = 2.2-4.5 eV (Athr - 
2240-3450 A). In the extremely metal-poor atmosphere, devia- 
tions of the mean intensity of ionizing ultraviolet radiation from 
the Planck function are more pronounced compared with those 
of the solar metallicity model (Fig.|8]l, resulting in much stronger 
departures from LTE. Figure|9]shows that all levels of Mg i, Ca i, 
and Fe i and the three lowest levels of Si i are strongly underpop- 
ulated in the line formation layers of the 5811/4.0/-4.5 model. 
Here we use the departure coefficients bi - n^^^^fnY^, where 
^NLTE ^jj^ ^LTE ^j-g j.jjg statistical equilibrium and thermal (Saha- 
Boltzmann) number densities, respectively. Non-LTE leads to a 
weakening of the Mg i. Si i, Ca i, and Fe i lines and positive non- 
LTE abundance corrections Anlte = logeNLXE - logfiLXE- We 
comment on the obtained results for individual species. 











i|f'ii!if.'iipf 



11 . l" ' 



2200 2400 

Wavelength [A] 



Fig. 8. Emerging fluxes calculated using the 5811/4/-4.5 (solid 
black) and 5780/4.44/0 (dotted violet) model. Bold marks indi- 
cate the thresholds of the Fe i levels with excitation energy of 0.9 
to 3.7 eV. 



The observed Mg i lines arise in the transitions 3/7 ^P° - 
3c/ ^D (382.9-383.8 nm) and 3/?^P° - 4s^S (517.2, 518.3 nm). 
For each line, the upper level is depleted to a lesser extent with 
regard to its LTE population than the lower level. Therefore, the 
line is weaker compared with its LTE strength not only because 
of the general overionisation {bi < 1), but also because the line 
source function {Si„ - bu/bi By) rises above the Planck function 
(By) in the line formation layers. Here, b,, and bi are the departure 
coefficients of the upper and lower levels, respectively. All inves- 
tigated lines have similar non-LTE abundance correction at the 
level of +0.2 dex from the calculations with 5h = 0. 1 ( Table|4]i. 
As expected, the departures from LTE are reduced for increased 
Hi collision rates {Su = !)■ 

The effect of bu/bi > 1 resulting in S u, > By is more 
prominent for the only available line of silicon. Si i 390.5 nm. 
Its lower level 3p'S follows the ground state of Sii inside 
logT5ooo < -1.5 due to collisional coupling, and it is strongly 
underpopulated in the line formation layers. Its coupling to the 
high-excitation levels for the upper level 4s 'P° turns out to be 
stronger than a coupling to the lower excitation levels, and tends 
to follow the continuum. Sin. This explains why Sii 390.5 nm 
has a larger non-LTE correction of Anlte = 0.34 dex (Sh = 0. 1) 
compared to the corresponding values for the Mg i lines and why 
Anlte is only slightly reduced when moving to S'e = 1 (Table|4]i. 
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Fig. 9. Depaiture coefficients h, for the selected levels of the in- 
vestigated atoms as a function of log rgooo in the model atmo- 
sphere 5811/4.0/-4.5. Bottom panel: we show every fifth of the 
first 60 levels for Fe i. They are quoted in the bottom-right part 
of the panel. All remaining higher levels of Fe i are plotted by 
dotted curves. We show the same for Fe 11. They are quoted in 
the top-right part of the panel. 5h = 0. 1 throughout. 



Fig. 10. Non-LTE abundance corrections for the Fe i lines in the 
5811/4.0/-4.5 model depending on fiexc from the calculations 
with 5h = 0.1 (filled circles) and 5h = 1 (triangles). 



For the resonance line of Cai at 422.6 nm, the non-LTE 
mechanisms are very similar to those of the Mg i lines. Calcium 
is the only element observed in SDSS J102915 + 172927 in two 
ionisation stages. Can dominates the element number density 
over atmospheric depths. Therefore, no process seems to affect 
the Can ground-state population, and 4j keeps its thermody- 
namic equilibrium value. The levels 3^/ and 4/7 follow the ground 
state in deep layers, and their coupling is lost at the depths out- 
side logTsooo < -1 where photon losses in the weakest line 
849.8 nm of the multiplet 3t/ - \p start to become important. 
In these atmospheric layers, bujbi < 1 is valid for each investi- 
gated line of Ca 11 resulting in dropping the line source function 
above the Planck function and enhanced line absorption. For the 
resonance line Can 393.3 nm, departures from LTE occur only 
in the very core and Anlte amounts to -0.07 dex. Non-LTE cor- 
rection is larger in absolute value for the IR lines of multiplet 
lid - 4p, 849.8, 854.2, and 866.2 because of the overpopulation 
of the lower level. 

Weakening of the Fe i lines is mainly due to overionisation. 
In SDSS J102915+172927, we measured only the low-excitation 
Fe I lines, with fiexc = 0-1.5 eV. The source function is quite sim- 
ilar to the Planck function for each investigated line, because all 
levels with Egxc = 0-4.5 eV behave similarly (Fig.|9]l- With the 
very similar behaviour of the departure coefficients for the lower 
levels, we calculated very similar non-LTE corrections, as can be 
seen in Fig. [10] Anlte varies between 0.29 and 0.36 dex in the 
calculations S'h = 0. 1 . Similarly to the Mg i lines, the departures 
from LTE are reduced significantly for S'h = 1 , see Table|4] 

Although only an upper limit was estimated for the Sr abun- 
dance, we performed the non-LTE calculations for Srii with 
[Sr/Fe] = -5.1. Non-LTE leads to a weakened Srii 407.7 nm 
line, and Anlte amounts to 0.16 dex if pure electronic collisions 
are taken into account in the SE calculations and decreases down 
to 0. 12 dex for 5h = L For Ti 11, we estimated a non-LTE correc- 
tion of -0.01 dex, assuming that the departures from LTE for the 
investigated Ti 11 lines are similar to those of the Fe 11 lines of 
similar excitation energy and equivalent width. 



5. The ISM towards the star SDSS J1 0291 5+1 72927 

The spectrum of SDSS J 1029 15 -1-172927 shows interstellar ab- 
sorptions of the Na i D-line doublet at 589.0 nm (see Fig. II lb and 
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Table 4. Non-LTE abundance corrections (dex) for the lines in SDSS Jl 029 15+ 172927 depending on surface gravity. 



Species 


A 


■C-exc 


Wob. 


5811/3.5/-4.5 


5811/4.0/-4.5 


5811/4.5/-4.5 




(nm) 


(eV) 


(pm) 


Sh = 0.1 


5h = 1 


5h = 0.1 


5h = 1 


5h = 0.1 


5h=1 


Mgi 


382.9355 


2.709 


0.8 


0.23 


0.08 


0.17 


0.04 


0.12 


0.02 


Mgi 


383.2300 


2.709 


W 


0.25 


0.09 


0.19 


0.05 


0.13 


0.03 


Mgi 


383.8290 


2.709 


bl 


0.28 


0.10 


0.21 


0.06 


0.14 


0.03 


Mgi 


517.2684 


2.712 


1.1 


0.23 


0.09 


0.19 


0.05 


0.13 


0.03 


Mgi 


518.3604 


2.717 


1.6 


0.24 


0.09 


0.19 


0.05 


0.13 


0.03 


Si I 


390.5523 


1.909 


1.8 


0.31 


0.30 


0.34 


0.30 


0.35 


0.27 


Cai 


422.6728 


0.0 


2.4 


0.32 




0.28 




0.24 




Can 


393.3663 


0.0 


27.7 


-0.08 




-0.07 




-0.07 




Can 


849.8023 


1.692 


2.2 


-0.12 




-0.11 




-0.09 




Can 


854.2091 


1.670 


8.6 


-0.36 




-0.30 




-0.22 




Can 


866.2141 


1.692 


7.0 


-0.28 




-0.24 




-0.18 




Fei 


385.9911 


0.0 


6.0 


0.40 


0.18 


0.32 


0.13 


0.23 


0.12 


Fei 


360.8859 


1.011 


2.8 


0.47 


0.19 


0.35 


0.13 


0.25 


0.13 


Fei 


407.1738 


1.608 


1.3 


0.43 


0.18 


0.33 


0.12 


0.24 


0.12 


Srn 


407.7709 


0.0 


none 


0.18- 


0.16 


0.16^ 


0.12 


0.142 


0.10 



' line is blended 

^Sh = 



the Caii-K line at 393.3 nm and Caii-H line at 396.8 nm (see de- 
tail inside Fig.|5]l. 

The interstellar feature is satisfactory modelled with one sin- 
gle component model providing a column density of log (Na i) 
= 12.11+0.01 cm-2 and log (Can) = 12.02 + 0.04 cm^^. The 
broadening of the lines is 7.3 ±1.1 kms"' in the Can lines and 
5.2 ± 0.1 kms^'in Nai, suggesting that the turbulence is the 
dominant broadening factor and that the two ions do not sample 
precisely the same material with the Can lines tracing ionised 
gas not detected in Na i. 

The Na i column density is consistent with that observed to- 
wards T] Leo, which at an angular distance of a few degrees 
shows log N(Nai)=12.08 cm^l 

High spectral resolution absorption observations of the lo- 
cal neutral and partially ionised interstellar gas using both the 
Na I D-line doublet and the Ca ii-K line allowed us to construct 
a picture of the distribution of neutral gas in the local interstel- 
lar medium. The Sun is placed within a low density interstellar 
cavity connected by interstellar tunnels to other surrounding cav- 
ities. 

Welsh et al. (l2010h reported on the results of an extended 
survey of Na i and Ca n absorption lines recorded at high spectral 
resolution towards thousands of early-type stars mostly located 
within 800 pc of the Sun. In particular, maps of the distribution 
of Na I absorption revealed that the local cavity has a 50 pc di- 
ameter and a 200 pc long extension in the direction of the star /3 
CMa and that there is an extension of rarefied gas into the lower 
galactic halo that forms an open-ended chimney feature. 

The line of sight towards SDSS J102915+172927 has the 
galactic coordinates 1-222.1 and b-56 i.e. lying in the third 
quadrant and at fairly high latitudes, intercepting low-density 
material. 

The observed column density ratio in the line of sight to- 
wards SDSS J102915+172927 is N(Nai)/N(Can)=l.l. This ra- 
tio is a diagnostic of the physical conditions in the diffuse in- 
terstellar gas. In cold (T~30K) and dense gas clouds most of 
the gas-phase Ca is depleted onto grains, and the Na i/Ca 11 ra- 
tio is >> 1. In warmer (T~1000K) and lower density ISM 
much of the Ca remains in the gas phase, and ratios of < 1.0 
are commonly found. For instance, for sight-line distances in- 




Fig. 11. IS lines of Nai in the spectrum of 
SDSSJ102915+172927. A telluric absorption Hne is de- 
noted by a crossed circle. The solid line is a synthesis of the IS 
one-component model. 



side the local cloud the Nai/Caii column density ratio for the 
warm (T~7000K) local interstellar gas clouds is ~ 0.2 . The 
Na i/Ca II ratio is thought to be also connected to the gas cloud 
velocity, such that at velocities > 30 km s ' interstellar dust 
grains may be destroyed by sh ocks and Ca is liberated into the 
gas phase (iRoutlv & Spitzer1[l952h . Accordingly, a low value 
of the Nai/Can ratio could be due either to the presence of 
warm and partially ionised gas and/or the presence of interstel- 
lar shocks. Since we do not observe high velocities, the ob- 
served low Na i/Ca 11 ratio is most probably caused by significant 
amounts of partially ionised Can gas. 

The Nai/Caii ratio value generally falls in the range 0.5 to 
20 but with a smaller range of ratio values, between 0.1 to 5, in 
quadrant 3 where our line of sight lies, and it matches the value 
measured towards SDSS J102915 + 172927. 

By using Nai as a tracer of neutral hydrogen, we can in- 
fer a total column density of log N(H) = 20.38 cm"' (Ferlet et 
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al. 119851) . This column density is similar to that of log N(H) - 
20.43 cm~^ directly measured in p Leo from the Lya absorp- 
tion U ne by means of lUE spectra (Shull & Van Steenberg 
Il985h . Because this is located at a distance of 950 pc, the gas 
column densities place a robust lower limit to the distance of 
SDSSJ102915+172927. 

With the given hydrogen column density the extinction is A^, 
=0.13, and assuming an extinction-to-reddening ratio of A„/E(B- 
V) = 3.1, we obtain E(B-V)= 0.04, which is quite consistent with 
the extinction value as deduced from extinction maps. 

6. Discussion 

The sample o f BMP candidates se lected from the SDSS archive 
(see details in lCaffau et al.l l201 lb), and observ ed with the high- 
efficien cy spectrographs X-Sho oter (.Caffau et alj|2011h) . or 
UVES dBonifacio et al.ll2011bl) is confirmed by this discovery 
to contain stars at the lowest levels of metallicity observed in 
the Galaxy. The few stars selected from SDSS and observed at 
high resolution cannot allow us to assess the real efficiency of 
the method, but this could be achieved with more observations. 

This discovery gives our understanding of the forma- 
tion of the first generations of stars a new perspective. The 
level of chemical enrichment of a star is characterised by 
the global metallicity Z. The most metal-poor stars formerly 
known had a metallicity of about Z - lO^^^Z©, while 
SDSSJ102915+172927 has a metallicity Z < IO-'^^Zq. A 
few stars with extrem ely low iro n abundances (of ffie order of 
[Fe/H1=-5.0 or lo wer, Christlieb et al. 2002; Frebel et al.ll2005l 
iNorris et al]|2007h were discovered in ffie past decade. However, 
their chemical composition is highly peculiar, with C, N, and O 
being more abundant than expected with respect to iron by sev- 
eral orders of magnitude, corresponding to a moderately low Z, 
of the order of lO^^Z©. This state of affairs supported the notion 
that the first generations of stars were formed exclusively by ex- 
tremely massive stars (now long extinct) and that there exists 
a critical metallicity be l ow which low-mass st ars cannot form 
dBronmi & Loebl l2003t ISchneider et~ai]|2003b . The existence 
of SDSS J 1029 15 -Hi 72927 demonstrates that either such a criti- 
cal metallicit v does not exist, in line wiffi the scenario of Greif 
et al. (2011), or that it is lower than 10"^Zq, in line with the 
scenario of Schneider et al. (2003). The peculiar chemical com- 
position of HE0107-5240 an d HE 1327-2326 is consistent with 
the pic ture by .Frebel et al.l (l2007h . using the iBromm & Loebl 
(l2003h ffieory, in which the critical parameter is not metallicity, 
but a suitable combination of C and O abundances. The fact that 
SDSS J102915-I-172927 is not carbon-enhanced places this star 
in what this theory calls the "forbidden zone". We have no upper 
limit on the oxygen abundance, but because C, N, and Mg are 
not enhanced (as they are in HE0107-5240 and HE 1327-2326), 
there is no reason to suspect any significant O enhancement. If 
we assume, conservatively, [O/H] < -4.1 and [C/H] < -3.8, the 
transition discriminanO is D < -3.7, while low-mass star forma- 
tion should only be possible for D > -3.5. Our discovery will 
a give new impulse to the formation theories of low-metallicity 
stars. The notion of critical metallicity or metallicity discrimi- 
nant may have to be revisited. 

The properties of SDSS J102915-H17 2927 are cornpatibl e 
with the protogalaxy scenario discussed in lOmukai et al.l (l2008l) . 
In this case, metal-poor gas under extreme UV irradiation, wiffi 
a normal dust-to-gas ratio and with a metallicity of Z> 5 x 
10"^ of Zq, fragments into a dense cluster of sub-solar mass 



stars (see ffieir figure 5). This particular way of forming low- 
mass metal-poor stars does not require any C/N/O enhance- 
ment, be cause o f the p eculiar thermal evolution of the gas. Also 
Schneid er et al.l (1201 ih expect formation of stars as metal-poor 
as SDSS J102915H-172927, and the cooling of the gas cloud is 
due to dust. In boffi cases a minimum metallicity is expected 
for low-mass star formation. On the other hand, the fragmen- 
tation found by the simulations of Clark et al. (2011) relies on 
H2 cooling and requires no metal enrichme nt. If the IMF of ffi e 
primordial stars is indeed flat, as found by iGreif et aP (1201 ih . 
ffien we may expect to find low-mass stars even more metal- 
poor than SDSS J102915+172927 and, possibly, of primordial 
metallicity. One concern on the observability of these primor- 
dial stars is the possibility that their atmospheres may have been 
polluted by metals during encounters with molecular clouds dur- 
ing their long lifetimes. While early estimations of this process 
predicted significant pollution to take place dYoshii et al.ll 19951) . 
subsequent observations fa iled to detect the predicted "accretion 
plateau" of light elements (iMolaro et al.lll997 ■ iBoesgaard et al.l 



D„ 



log(l0lc/HJ + 0.3 X IQlo™) 



1999';'Smilianic et al."2009';'Duncan et ari997'; Garcia Lopez et 
al. 1998). More recent estimations of the accretion process effi- 
ciency (Frebel et al. 2009) seem to exclude any significant pol- 
lution, es£ecmn^MfJhe_gmnordial stars have a weak solar-like 
wind d Johnson & Khochfaill201ll) . These considerations should 
encourage an extensive search for more stars that are as metal- 
poor as SDSS J102915+172927 or even more metal-poor. 

Increasing the number of stars at these very low metallici- 
ties will also help understanding the behaviour of lithium at ex- 
tremely low metallicities. If below a given metallicity low-mass 
stars can only be formed by fragmentation of larger collapsing 
clouds (Clark et al. 2011; Greif et al. 2011), a distinct possibil- 
ity is that this process results in turbulence and/or rotation fast 
enough to destroy the primordial lithium originally present in 
the gas effectively during the pre-main sequence phase. In this a 
scenario we expect all stars below a given metallicity to show no 
lithium. As star formation begins to proceed in a more conven- 
tional way, through metal-line cooling, the stars do not destroy 
the lithiu m, and the Spite plate au appears. The Spite plateau 
meltdown dSbordone et al.ll2010h could simply mark a transition 
region between ffie two regimes of star-formation. 
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Appendix A: Variability (absence of) 

To verify possible magnitude variations of 
SDSSJ102915-I-172927 either on a long or a short time 
scale we obtained a sample of images at the Asiago Schmidt 
Telescope and also checked if other images in public databases 
were available. In the last case we only found four images in 
the Isaac Newton Telescope (INT) database as part of the Wide 
Field Survey (WFSur) of the Northern Hemisphere. 



A.1. Asiago Schmidt Telescope 

On the night of 5 May 201 1 we obtained a sequence of 4x300sec 
images in the V and 4 x 300sec images in the Cousins-R band 
at the 92/67 Schmidt Telescope at Cima Ekar Observatory in 
Asiago equipped with a KAI-llOOOM 4.0 x 2.6 milHon pixels 
CCD. With a pixel size of 0.85" the camera has a field of view of 
56x30 arcmin^. The seeing was not particularly good, ^1.9 arc- 
sec, but it allowed us to obtain data with a good sampling of the 
Point Spread Function. The data also served as a second epoch 
for the proper motion measurements (see below). We corrected 
the raw images in the usual way with appropriate master bias 
and flat-fields. We then used DAOPHOT (Stetson 1987) for ex- 
tracting PSF photometry from the images. The calibration of the 
data was performed in a relative way matching all the sources in 
the field with the SDSS photometry of the star The R band used 
here is not too far from the SDSS r filter, therefore we calibrated 
it against the r, {r - i) SDSS magnitude and colour We obtained 
a mild colour term of -0.09 + 0.01. To calibrate the V band, 
which is not available from SDSS, we used a linear combination 



21 



Caffau et al.: An extremely primordial star 



of the SDSS g and r band appropriate for the transformation to 
the Johnson-Cousin system for main sequence stars, which can 
be found on the SDSS web pagqj: the magnitude calculated in 
this way is V= 16.688 + 0.005. 

The results are reported in Table lA.fl and shown in Fig. lA.ll 
where in the top panel the four V band measurements and in the 
bottom the four R band magnitudes are shown. As a reference 
the SDSS magnitude of SDSS J102915+172927 is shown as re- 
ported in Table [1] The average errors of the Cima Ekar Schmidt 
measurements are ^ 0.035 mag for the R band and ^ 0.032 mag 
for the V band. The systematic difference from the SDSS pho- 
tometi-y is r= -1-0.006 mag and for V= -1-0.004. 

A.2. INT La Palma 

Only four images in the r band are present in the INT archive 
of images for the Wide Field camera (WFC). These images are 
part of the WFSufl, a collection of several programmes aim- 
ing at providing a quick coverage of the northern sky either in 
single- or multi-band. The PI of the observing programme is N. 
Walton, and to our knowledge the final catalogue and images 
have not yet been released but only intermediate products. The 
area surrounding our star was imaged on the night of 10 March 
2003 with four images in the r band of 1 minute each. We bias- 
and flat-field corrected the images using material from the same 
night and produced a DAOPHOT photometry of only the CCD 
where the star is present (CCD nr 3). We then calibrated the raw 
WFC photometry using the SDSS photometry of all stars in the 
field, properly excluding saturated and faint stars. The results 
are shown in the right column of Fig. lA.ll in which we show the 
variation of the magnitude in time. The plot shows the reference 
SDSS r-band magnitude as listed in TablT] The scatter of the 
measurements is ^ 0.007 mag with a 0.000 mag difference with 
the SDSS value. 

With the material at our disposal we can safely conclude that 
there is an absence of long-term photometric variations at the 
level of Ar < 0.006 variations considering the errors of the star. 
For the short term the variations seems to be insignificant. A 
longer-time monitoring with a better timing is underway to con- 
firm this overall trend, considering that the monitoring covered 
less than 1 hour of data-points for the Asiago Schmidt and 12 
minutes for INT data. 



Appendix B: Atomic data line-by-line 
References 

Stetson, P. B. 1987, PASP, 99, 191 
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Fig.A.l. Magnitude measurements for SDSS J102915-H172927 
from Asiago and La Palma. On the left column are shown the 
Asiago Schmidt telescope measurements versus time (in hours) 
for the night of 5 May 201 1 . Top panel, the 4 V magnitude com- 
pared to a linear combination of the SDSS g and r magnitudes. 
Bottom panel, the 4 R magnitude measurements calibrated vs the 
SDSS r band. In both panels the respective SDSS magnitudes are 
shown as short dashed lines and Icr errors (long-dashed lines). 
On the right column are shown the r band measurements of the 
WFC/INT on La Palma for the night of 10 March 2003. 

Table A.l. r and V photometry from Asiago Schmidt and from 
INT for SDSS J102915+172927. 



Asiago/Schmidt 


MJD 


Band 


Magnitude 


55686.919175 


r 


16.523 ± 0.035 


55686.927659 


r 


16.544 ± 0.037 


55686.935613 


r 


16.549 ± 0.035 


55686.943686 


r 


16.528 ±0.038 


55686.923177 


V 


16.684 ±0.032 


55686.931663 


V 


16.672 ±0.032 


55686.939619 


V 


16.702 ±0.031 


55686.947691 


V 


16.676 ±0.033 




WFC/INT 


MJD 


Band 


Magnitude 


52708.897580 


r 


16.548 ± 0.003 


52708.898767 


r 


16.547 ± 0.002 


52708.899959 


r 


16.533 ± 0.002 


52708.901151 


r 


16.542 ± 0.002 



* http://www.sdss.org/DR7/algorithms/sdssUBVRITransform.html 
^ http://www.ast.cam.ac.uk/ wfcsur/index.php 
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Table B.l. 


Lines analysed in this work. 






Element 


A 


tl'low 


log 9f 


[X/H]lD_L7£ 


Notes 


Mgi 


382.9355 


2.709 


-0.207 


-4.869 




Mgi 


383.2 


4.346/2.712 


-0.456/+0.021 


-4.574 


blend 


Mgi 


383.8290 


2.717 


-1.583/+0.339/-0.408 


-4.760 


blend 


Mgi 


517.2684 


2.712 


-0.402 


-4.669 




Mgi 


518.3604 


2.717 


-0.158 


-4.699 




Si I 


390.5523 


1.909 


-1.090 


-4.270 




Cai 


422.6728 


0.000 


+0.243 


-4.718 




Can 


849.8023 


1.692 


-1.416 


-4.719 




Can 


854.2091 


1.670 


-0.463 


-4.500 




Can 


866.2141 


1.692 


-0.723 


-4.593 




Tin 


334.9408 


0.049 


+0.430 


-4.819 




Tin 


336.1218 


0.028 


+0.280 


-4.795 




Tin 


337.2800 


0.012 


+0.180 


-4.432 




Tin 


338.3768 


0.000 


+0.142 


-4.723 




Tin 


375.9296 


0.607 


+0.270 


-4.791 




Tin 


376.1323 


0.574 


+0.170 


-4.964 




Nil 


339.2983 


0.025 


-0.540 


-4.363 




Nil 


341.4760 


0.025 


-0.060 


-4.422 




Nil 


344.6255 


0.109 


-0.410 


-4.558 




Nil 


345.8456 


0.212 


-0.260 


-4.662 




Nil 


346.1649 


0.025 


-0.360 


-4.367 




Nil 


349.2954 


0.109 


-0.270 


-4.671 




Nil 


351.0332 


0.212 


-0.670 


-4.627 




Nil 


351.5049 


0.109 


-0.260 


-4.670 




Nil 


352.4535 


0.025 


-0.030 


-4.416 




Nil 


361.9386 


0.423 


-0.040 


-4.710 




Fei 


355.8515 


0.990 


-0.629 


-4.846 




Fei 


356.5379 


0.958 


-0.190 


-4.846 




Fei 


357.0098 


0.915 


+0.153 


-4.932 




Fei 


358.1193 


0.859 


+0.406 


-4.780 




Fei 


360.8859 


1.011 


-0.100 


-4.842 




Fei 


361.8768 


0.990 


+0.000 


-4.847 




Fei 


370.5566 


0.052 


-1.334 


-4.522 




Fei 


370.9246 


0.915 


-0.646 


-4.595 




Fei 


371.9935 


0.000 


-0.431 


-4.466 




Fei 


372.7619 


0.958 


-0.631 


-4.642 




Fei 


373.7132 


0.052 


-0.574 


-4.645 




Fei 


374.5561 


0.087 


-0.771 


-4.633 




Fei 


374.5899 


0.121 


-1.335 


-4.660 




Fei 


374.8262 


0.110 


-1.016 


-4.720 




Fei 


375.8233 


0.958 


-0.030 


-4.799 




Fei 


376.3789 


0.990 


-0.240 


-4.725 




Fei 


376.7192 


1.011 


-0.390 


-4.812 




Fei 


378.7880 


1.011 


-0.860 


-4.702 




Fei 


381.5840 


1.485 


+0.240 


-4.814 




Fei 


382.0425 


0.859 


+0.120 


-4.819 




Fei 


382.4444 


0.000 


-1.360 


-4.581 




Fei 


382.5881 


0.915 


-0.040 


-4.901 




Fei 


382.7822 


1.557 


+0.060 


-4.954 




Fei 


383.4222 


0.958 


-0.302 


-4.782 




Fei 


384.0437 


0.990 


-0.506 


-4.778 




Fei 


384.1048 


1.698 


-0.050 


-4.771 




Fei 


384.9966 


1.011 


-0.970 


-4.590 




Fei 


385.6371 


0.052 


-1.290 


-4.647 




Fei 


385.9911 


0.000 


-0.710 


-4.592 




Fei 


388.6282 


0.052 


-1.080 


-4.591 




Fei 


389.5656 


0.110 


-1.670 


-4.611 




Fei 


389.9707 


0.087 


-1.530 


-4.579 




Fei 


392.2912 


0.052 


-1.650 


-4.648 




Fei 


392.7920 


0.110 


-1.590 


-4.602 




Fei 


393.0297 


0.087 


-1.590 


-4.463 




Fei 


404.5812 


1.485 


+0.280 


-4.831 




Fei 


406.3594 


1.557 


+0.070 


-4.851 




Fei 


407.1738 


1.608 


-0.020 


-4.823 




Fei 


427.1760 


1.485 


-0.160 


-4.830 




Fei 


430.7902 


1.557 


-0.070 


-4.899 




Fei 


432.5762 


1.608 


-0.010 


-4.884 




Fei 


438.3545 


1.485 


+0.200 


-4.735 




Fei 


440.4750 


1.557 


-1.400 


-4.760 




Srii 


407.7709 


0.000 


+0.167 


-5.139 
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